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Abstract 

We report on a detailed investigation of the next-to-leading order (NLO) QCD 
corrections to H^7+jet production at the Tevatron and the LHC using a fully-flexible 
parton-level Monte Carlo program. We include the full leptonic decay of the W, tak- 
ing into account all off-shell and finite width effects, as well as non-standard WWj 
couplings. We find particularly sizable corrections for the currently allowed parameter 
range of anomalous couplings imposed by LEP data. In total the NLO differential 
distributions reveal a substantial phase space dependence of the corrections, leaving 
considerable sensitivity to anomalous couplings beyond scale uncertainty at large mo- 
mentum transfers in the anomalous vertex. 
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1 Introduction 



Electroweak diboson production in association with a hard jet is an important class of 
processes at hadron colliders such as the Large Hadron Collider (LHC) or the Tevatron. 
Electroweak boson phenomenology generically provides a window to the electroweak sym- 
metry breaking sector, and diboson signatures with or without jets are therefore potentially 
sensitive to new interactions beyond the Standard Model (BSM). Equally important, SM- 
diboson+jet production contributes to the irreducible background of new physics searches 
in these channels. Hence, precise cross section predictions are mandatory to obtain a correct 
interpretation of possible excesses, which might be observed in the near future. The total 
cross sections for the diboson+jet production processes are fairly large compared to the pure 
diboson production channels if extra jet emission at large available center-of-mass energy is 
kinematically unsuppressed. The large one jet-inclusive cross section is mainly due to ac- 
cessing the (anti)proton's gluon parton distribution function at small momentum fractions 
already at leading order (LO). At the same time, the gluon- induced partonic subprocesses 
opening up at 0{a s ) give rise to exceptionally large next-to-leading order (NLO) QCD cor- 
rections to inclusive diboson production, see Ref. [1]. Qualitatively similar observations have 
also been made for the NLO QCD corrections to various diboson+jet production processes 
in a series of recent publications [2-9]. 

In the present paper we extend our NLO calculation of pp,pp — > W ±/ y + jet + X [5] to 
anomalous WW^/ couplings*. We review the SM W ±r y + jet phenomenology in detail and 
discuss its modifications due to anomalous couplings at NLO QCD precision. This allows 
us to discriminate between the effects of new physics in terms of effective interactions from 
the impact of higher order corrections. Anomalous couplings searches represent benchmark 
tests for non-SM interactions at the LHC at small integrated luminosity C < 30 fb _1 (see, 
e.g., Ref. [10]). Measurement and discovery strategies have received lots of attention, both 
from the theoretical (e.g. Refs. [11-14]) and the experimental side (e.g. Refs. [10,15-18]). 
In this context, diboson production processes are important channels at the LHC because 
they exhibit large total rates, and, in case of W ±r y production, because they are sensitive to 
deviations of the underlying electroweak model from the SM via so-called radiation zeros. 
These classical zeros of the amplitude in the qQ — > W ±r ) channels at the photon center-of- 
mass scattering angles cos# 7 = +1/3 are special to the completely destructive interference of 
gauge boson-radiation in an unbroken renormalizable field theory, Ref. [19]. Any deviation 
from QED by additional non-SM operators ultimately destroys this characteristic radiation 
pattern. At the LHC, the antiquark direction is, in principle, indistinguishable from the 
quark direction because of the proton-proton initial state, and the radiation zero gets con- 
siderably washed out. "Signing" the quark direction according to the event's overall boost, 
which has been considered in the context of dilepton asymmetries and electroweak mixing 
angle measurements [20], has been shown to efficiently lift the initial state's degeneracy in 
Ref. [10]. Furthermore, the radiation zero remains present only if additional electromagnet- 
ically neutral (e.g. gluonic) radiation is collinear to the photon. Hence, additional QCD 

*For convenience we refer to the computed processes as W ± j + jet production even though we include 
all finite width and off-shell effects of the massive W. 
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emission, as part of the NLO contribution to W'y production, is dangerous to observing the 
radiation zero. At the same time, the radiation zero becomes nearly impossible to measure 
in W^7+jet production [13]. 

Crucial to significance-improving strategies [11] is therefore an additionally-imposed jet 
vetoL Jet vetoing, however, is a delicate strategy in fixed-order perturbation theory from a 
theoretical point of view. The observed reduction of scale dependence for the exclusive Wy 
production at the LHC (see Ref. [1]) is predominantly due to excluding a region of phase 
space from the total inclusive cross section, which is well-accessible at the large available 
center-of-mass energy. This is also reflected in additional jet radiation becoming highly 
probable as part of the real emission contribution to the NLO diboson cross section, a(W , y + 
jet)/cr(H^7) ~ 3. Hence, the dominant perturbative uncertainties result from the W'y + jet 
contribution, which is a leading order a s contribution to NLO Wy production. Consequently, 
current Monte Carlo-driven strategies that involve jet vetos to measure anomalous couplings 
from fits to high transverse momentum distributions (via e.g. neuronal net algorithms trained 
to the NLO Wy distributions) inherit significant uncertainties, considerably larger than 
those given by scale variations of the exclusive NLO cross sections. By computing Wy + jet 
production at NLO accuracy, we are able to realistically estimate the anomalous parameters' 
impact on the vetoed cross section and contribute a crucial part towards modelling inclusive 
Wy production at a higher perturbative precision. 

We organize this paper as follows: Section 2 gives details on our Monte Carlo implemen- 
tation and introduces the notion of anomalous WWy couplings to the reader. In Sec. 3, we 
discuss the numerical results; we give total cross sections and differential distributions, both 
for SM and anomalous production. We also comment on the sensitivity to anomalous cou- 
plings in NLO QCD W'y + jet production at the LHC and we quote W'y + jet cross sections 
for selection cuts adapted to anomalous couplings' searches. The phenomenological impact 
of anomalous couplings on H^7+jet production at the Tevatron is too small for parameter 
choices that are compatible with the bounds imposed by LEP data, Ref. [17]. Hence, we only 
quote Tevatron results for SM-like production in Sec. 3.3. Section 4 closes with a summary 
and gives an outlook to future work. Our Monte Carlo code will become publicly available 
with an upcoming update of Vbfnlo [21]. 



2 Details of the calculation 



There are three contributing partonic subprocesses at 0(a 3 a s ) for pp,pp — > e v e yj + X, 

qQ — li'vag, (la) 
Qg — >Cv l yq, 

qg — > fvaQ, (lc) 



^The jet veto also removes kinematical configurations which are less sensitive to anomalous couplings 
due to small momentum transfers in the WWy vertex from the total cross section, see below. 
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Figure 1: Representative Feynman graph 
contributing to the virtual corrections to 
the partonic subprocess ud — > e~u e jg at 
0(a 3 a 2 ). The crosses mark points where 
the photon can be attached to the quarks 
or the W for the fixed gluon-PF order along 
the quark line. 



not counting the subprocesses, which follow from interchanging the beam directions. We use 
the shorthand notation q = (d, s) and Q = (u, c) and assume a diagonal CKM matrix. At 
the LHC, a non-diagonal CKM matrix decreases our leading order results only at the per 
mil-level. Unitarity of the (non-diagonal) CKM matrix guarantees that all CKM-dependence 
drops out for flavor-blind observables computed from the dominant gluon-induced subpro- 
cesses. At the Tevatron, we find our cross sections decreased by about 3%. Both modifica- 
tions are smaller than the residual scale dependence at NLO, so that a diagonal CKM matrix 
is an adequate approximation for our purposes. Bottom quark contributions are absent at 
LO for the above approximations and can be further suppressed experimentally by b vetoing, 
and we therefore neglect bottom contributions throughout the computation. 

The LO matrix elements of Eq. (1) are calculated using Helas routines [22] generated 
with MadGraph [23]. Although we refer to the processes as H /± 7+jet production for 
convenience, we include all off-shell and finite width effects of the PF's decay to leptons, as 
well as the photon's coupling to the charged final state lepton and calculate the full QCD 
corrections to the processes pp,pp — > t~v^j + X and pp,pp — > t^v^j + X at 0{a z a 2 s ). 
Representative virtual Feynman graph topologies are sketched in Fig. 1. We perform the 
numerical phase space integration with a modified version of Vegas [24], which is part 
of the Vbfnlo package. We divided up the integration and explicitly sum over different 
channels that are optimized for the two and three-body decay of the (off-shell) W boson, 
PPiPP — > W^j —> £isjj and pp,pp — > Wj — > ivjj, respectively. For additional details on 
the phase space integration's validation against MadEvent [23] and Sherpa [25], we refer 
the reader to our previous publication [5]. The numerical implementation includes the finite 
width of the W by using the fixed width scheme* of Ref. [26]: The weak mixing angle is 
taken to be real and we use Breit-Wigner propagators for the massive W throughout. 

The counter term-renormalized virtual amplitude for, e.g., ud — > l~v e ^g (Fig. 1) in 
*This is also the width scheme used by MadGraph. 
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conventional dimensional regular izat ion d = 4 — 2e can be cast into the form, 




with Cf = 4/3 and Ca = 3 denoting the casimirs of the fundamental and adjoint represen- 
tations, respectively. The virtual amplitude exhibits an identical structure in color space as 
the Born matrix element and we implicitly assume the su(3) generator in the fundamental 
representation to be part of the definitions of .Mlo and -Mvirt- The constants 7$ are fixed 
by the SU (3) representations of the Np = 5 active quark flavors and the gluons, 

% = l~ q = \C F , % = ^C A - ^T r N F , (3) 

where T r — 1/2 is the Dynkin index of the fundamental representation. s,t,u denote the 
familiar Mandelstam variables of a 2 — > 2 process, taken to be space-like. Depending on the 
subprocess, the analytical continuation to physical kinematics is performed automatically 
within our numerical code by effectively restoring the propagators' e description. To arrive 
at the correct logarithms when evaluating Eq. (2) for physical s,t,u, we have to write, e.g., 

log(-s)=log(|s|)-i7re(-s) (4) 

for the qQ induced subprocesses for which s is a time-like quantity. Similar formulae have to 
be taken into account for the analytical continuation of the dilogarithms and can be inferred 
from the literature, e.g. from Ref. [27], appendix C. A^virt in Eq- (2) represents finite con- 
tributions that embrace tensor coefficients and fermion chains after algebraic manipulations. 
We also implicitly include the so-called rational terms to our definition of .Mvirt- These 
finite contributions arise from the interplay of (/-dimensional numerator and denominator 
algebra in the limit e — > (see, e.g., [28]). To compute the tensor coefficients, we apply the 
Passarino-Veltman recursion up to box topologies [29] and the Denner-Dittmaier reduction 
for pentagon graphs [30]. 

The full virtual amplitude can be assembled from elementary building blocks, that divide 
the one loop graphs into certain groups. This strategy has already been applied to calculate 
a series of different processes at NLO QCD precision in Ref. [31], both in the SM and beyond. 
Concretely, we sum all self-energy, triangle, box and pentagon corrections to a quark line 
with three attached gauge bosons of a given order to yield a single numerical routine. A 
similar routine is constructed from a quark line with two attached gauge bosons (see Fig. 1). 
These routines are set up with an in-house framework that partly uses FeynArts [32] 
and FeynCalc [33]. From these building blocks we construct the full SM loop amplitude 
for a given subprocess by trivial permutations of the bosons' momenta and polarization 
vectors, which also encode the two and three-body decay of the massive W depending on 
the building block. Replacing the SM polarization vectors by polarization vectors modified by 
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Figure 2: Sample Feynman graph con- 
tributing to the partonic real emission sub- 
process ud — > e~v e ^gg at 0(a 3 a 2 ). The 
gluon can also be attached to the quark and 
gluon lines at the positions marked by the 
circles. Feynman graph topologies, where 
the photon is radiated off at different posi- 
tions analogous to Fig. 1, are not shown. 



including anomalous WW'y couplings, we straightforwardly generalize our NLO calculation 
to anomalous IF7 + jet production. The renormalization in Eq. (2) is performed on-shell for 
the wave functions and we use the MS-scheme to renormalize the strong coupling constant. 
The virtual corrections to all other subprocesses of Eq. (1) can be recovered from Eq. (2) by 
crossing and analytical continuation analogous to Eq. (4). 

To speed up the numerical evaluation of the 90 subprocesses of the real emission contri- 
bution, we compute the respective matrix elements using optimized code that employs the 
spinor helicity formalism of Ref . [34] . The generalization to matrix elements with anomalous 
couplings is again performed by modifying the three-body effective W polarization vector to 
include the non-SM interactions of Sec. 2.1. The subprocesses for pp,pp — > e~v e r yjj + X at 
0(a 3 a 2 ) can be classified, modulo crossing, flavor summation, and initial state interchange, 
into 

mm — > t~vf^du , uu — > t~vf^sc , (5a) 
dd — > t~vf>fud , dd — > i~vf^cs , (5b) 
ud — ffvagg. (5c) 

Representative Feynman topologies for the last line's subprocess are indicated in Fig. 2. We 
store intermediate numerical results common to all subprocesses and reuse them whenever 
possible to speed up the numerical implementation. All matrix elements have been checked 
explicitly against code generated with MadGraph as well as against Sherpa for integrated 
cross sections (cf. [5]). Applying the dipole subtraction of Catani and Seymour [35], we have 
verified our implementation against code generated with the MadDipoles [36] package. Our 
code is optimized such that intermediate dipole results are reused in order to avoid redundant 
dipole or Born-level matrix element calculations. We also recycle the dipoles' Born-level 
matrix elements into the integration of the finite collinear remainder, which is left after 
renormalization of the parton distribution functions [35,38]. We integrate this contribution 
over the real emission phase space applying the phase space mappings of Ref. [39]. The 
remaining IR- singularities of the virtual matrix element Eq. (2) cancel analytically against 
the one-parton phase space-integrated dipoles, symbolically denoted by (J) in the language 




5 



of Ref. [35]. For, e.g., the gQ-induced channels adding the (I) operator yields 
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where we have already performed the analytical continuation to time-like s. The primed 
A^virt i R Eq. (6) indicates that we deal with a changed finite piece compared to Eq. (2) due 
to the analytic continuation. Note that we again include the Born level color structure into 
the definitions of the amplitudes; summing and averaging over colors and spins yields trivial 
subprocess-dependent additional prefactors. Equation (6) is explicitly free of IR divergencies 
and hence finite for e — > 0. A similar factorization and cancellation up to box diagrams has 
been demonstrated in Ref. [40]. 



2.1 Anomalous WW7 couplings 

We parametrize deviations of the SM electroweak sector by extending the SM Lagrangian 
by the most general, QED-preserving, Lorentz and CP-invariant operators up to dimension 
six [41], 



'WW-y 



— ie 



m w 



(8) 



It is customary to express the /t-induced deviation from the SM operators 5 by An = k — I. 
The parameters Ak and A are related to the electric quadrupole moment Qw and magnetic 
dipole moment /iw of the W boson [41], 



Hw = 77- — (2 + Ak + A), Q w = — ^-(1 + Ak-A). (9) 

2m w m A w 



Retaining unitarity at high energies is crucial to meaningfully modeling physics beyond 
the SM. On the one hand, if probability conservation is violated, the cross section receives a 
sizable contribution from probing the matrix elements at large invariant masses, even though 
the parton luminosities tame the matrix elements' unphysical growth in this particular phase 



§We recover the electroweak part of the SM by choosing the parameters k = 1 and A = 0. 
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space region. On the other hand, if unitarity is conserved, the phenomenology is mainly 
dominated by comparatively low invariant masses by the same reason. In order not to violate 
unitarity, the parameters Ak, and A have to be understood as low-energy form factors, and 
their precise momentum dependence does depend sensitively on physics beyond the SM. 
However, a widely-used phenomeno logical parametrization is (cf. Refs. [11,41]), 

Ak= {i+^+pwyi^y^ A= (i+(p 7 +^) 2 M 2 r' (10) 

where A represents the scale at which the beyond-the-SM interactions become strong, i.e. 
the scale at which W compositeness is resolved. p 1 and pw denote the final state momenta of 
the photon and the W, respectively. Unitarity imposes n K > 1/2 and ri\ > 1 [42]; customary 
choices by experimentalists are dipole profiles n K = ri\ = 2 [10,15-18]. Note that we do not 
include anomalous CP-violating operators since they have already been tightly constrained 
by measurements of the neutron electric dipole moment [11,43]. 

We include the anomalous interactions by constructing purpose-built Helas routines that 
implement the effective Lagrangian of Eq. (8) in a straightforward way using FeynRules 
(Ref. [44]) to arrive at the analytical expression for the anomalous vertex function. 



3 Numerical results 



3.1 Selection criteria, general Monte Carlo input and photon iso- 
lation 

Throughout, we use CTEQ6M parton distributions [45] with a s (mz) = 0.118 at NLO, 
and the CTEQ6L1 set at LO. We choose m z = 91.188 GeV, m w = 80.419 GeV and 
G F = 1.16639 x 10~ 5 GeV -2 as electroweak input parameters and derive the electromagnetic 
coupling a and the weak mixing angle from SM-tree level relations. The center-of-mass 
energy is fixed to 14 TeV for LHC and to 1.96 TeV for Tevatron collisions. We consider one 
family of light leptons in the final state, which we treat as massless, i.e. we quote results for 
pp,pp — > e~V e ~fj + X or pp,pp — > ji~Ji e ^i + X when we speak of W™7 + jet production. 

Jets are recombined from partons with pseudorapidity |^| < 5 applying the algorithm of 
Ref. [46] with resolution parameter D = 0.7. The reconstructed jets are required to lie in 
the rapidity range 

|%|<4.5. (11) 

The charged lepton and the photon are required to fall into the rapidity coverage of the 
electromagnetic calorimeter, i.e. we impose 

M<2.5, |r? 7 |<2.5. (12) 



7 



t- 

of 
1 



200 



^ 150 



100 



50 




0.2 0.4 0.6 0.8 1 1.2 1.4 



min i? 



800 




200 



r-, 150 



b 

X3 



100 



50 




0.2 0.4 0.6 0.8 1 1.2 1.4 



min R 




min Ra~ 



Figure 3: Minimal photon-jet separation in the azimuthal angle - pseudorapidity plane for 
W~j+iet production at the LHC. Plotted are the leading order and next-to-leading order dis- 
tributions for different isolation scales £ = pj,, p^'" 1 , max(pj,,^,^ ) in Eq. (13). The upper row 
displays the lower row's min ii,' 7 distributions around the isolation cone in more detail. The iso- 



lation parameter is 5q = 0.6 (indicated by the vertical line), and \xp 
chosen as described in the text. 



Hr = 100 GeV. Cuts are 



In order not to spoil the cancellation of the IR singularities and to minimize contributions 
from non-perturbative jet fragmentation associated with collinear photon-jet configurations, 
we apply the photon isolation criterion of Ref . [47] : 



partem, i 

Pt — 



1 — cos R 
1 — cos 5q 



(13) 



The index i in Eq. (13) runs over all partons found in a cone around the photon of size R 
in the azimuthal angle-pseudorapidity plane. The IR-safe cone size around the photon is 
given by 5$, and, in principle, can be replaced by an arbitrary energy scale E , which then 
determines the penetrability of the photon cone by soft QCD radiation. 



In a realistic experimental setting, the photon identification efficiency depends on the 



8 



number of finite-sized calorimeter towers that enter the candidate photon reconstruction (cf. 
Ref. [48]). It is therefore worthwhile to note that the isolation scale choice implicitly enters 
in the photon definition and different scale choices will be accompanied by different identi- 
fication efficiencies for small values of 5 . A detailed investigation of this direction should 
include effects ranging from pile-up and underlying event to NLO-corrected jet fragmentation 
and is beyond the scope of this work. For values of Sq much larger than the electromagnetic 
calorimeter cell size of AR ps 0.04, however, Eq. (13) is a "sliding-cut" prescription, which is 
experimentally well-defined on the level of already reconstructed particles. We can therefore 
compare the impact of the IR-safe isolation criterion to the NLO scale uncertainty, which 
turns out to be of order 10%. Replacing p^ in Eq. (13) by other intrinsic scales to the 
process, e.g. by maxp J T , corrects our NLO results at the level of 1% for inclusive cuts and 
5q = 0.6 (Fig. 3). Comparing the LO and NLO inclusive distributions we find a large net 
increase of the differential cross section. The differences in the choice of the isolation scale, 
however, are only visible in the jet-photon separation distribution around the photon cone 
mm Rj y PS 5 . For larger separations mmi?j 7 > S we do not find any notable phenomeno- 
logical impact of the isolation scale once we take into account the cross section's residual 
scale dependence on //R,//p. For the isolation scale S = p^T' 7 = \pj, + p? + p T \ the threshold 
behavior around 5q changes most significantly when comparing LO and NLO distributions 
(we denote the neutrino's four- momentum by i> in the following). The behaviour at LO is a 

consequence of p^' 1 = pip and the W recoiling against the photon-jet pair if the jet is emitted 
around the photon cone. Therefore, probing smaller parton-photon separations effectively 
means increasing the transverse momentum of the W for central events. The isolation scale, 
however, is set by the jet itself. Due to the exponential drop-off of the p^ = \p T + ft T \ spec- 
trum, collinear jet-photon configurations are highly attenuated for separations smaller than 
So at LO, see Fig. 3. At NLO the kinematical LO correlation of the jet and W^-system is 
modified by additional parton emission, which allows the W to be emitted at smaller trans- 
verse momenta. Thus, QCD radiation into the photon cone around the threshold S becomes 
more probable than at LO. At the same time, however, p^' 7 decreases and more partons get 
vetoed at distances smaller than So, and a steep drop-off is still visible in mmRj^ around 5q 
at NLO for S = pJ^' 7 in Fig. 3. In addition, we add events with a positive-definite weight to 
the minimum-separation distribution with the second resolved jet approximately balancing 
the jet-photon- W system in p T , which also modifies the mini?^ threshold behaviour. 

We now move on to investigate the general features of VF7+jet production with inclusive 
cuts on jets, photon, and lepton; we require 

p> T > 50 GeV, p e T > 20 GeV, p T > 20 GeV . (14) 

To avoid the collinear photon-lepton configurations we impose a finite separation in the 
azimuthal angle-pseudorapidity plane of 

i2, 7 = ^A$ 7 + A77f 7 >0.2. (15) 
For the jet-lepton separation we choose 

R jt > 0.2, (16) 
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and for the photon isolation we impose 



<5 = 0.6. (17) 

It is customary to also analyze the cross sections' behavior with an additional 'no resolvable 
2 nd jet '-criterion [2], i.e. a veto on the second jet, if it gets resolved, 

y/ cto > 50 GeV, \y] cto \ < 4.5 (exclusive NLO) . (18) 

This allows us to identify the dominant contributions to the NLO-inclusive cross section. In 
addition, it was shown in Ref. [5] that this cut leads to seemingly stable exclusive W^y+jet 
production cross sections at the LHC. Similar observations have been made for the other 
diboson+jet cross sections provided in Refs. [2-4,6-8]. The stabilization of the exclusive 
cross section does, however, not provide a reliable estimate of the perturbative cross sections' 
uncertainties over the whole phase space. We will discuss this in more detail in Sec. 3.2. 



3.2 NLO QCD W7+jet production in the SM at the LHC 

We quote the inclusive cross sections at the LHC for the chosen cuts in Tab. 1. The dif- 
ferences comparing W + to W~ is mainly due to the different parton distributions at the 
LHC. The cross sections' scale dependencies, estimated by varying /ip = (J,r by a factor 
two around 100 GeV, is reduced from approximately 11% to 7% when including the NLO 
QCD corrections. As already pointed out, introducing a veto on the second resolved jet, 
leads to considerable stabilization of the integrated NLO cross section; the one jet-exclusive 
cross sections (Tab. 2) exhibit scale variations of 0.5% for Vy - 7+jet and 2% for iy + 7+jet 
production, respectively. Naively, this suggests that vetoing the second hard jet amounts 
to a perturbative stabilization of the NLO cross section at a lower exclusive rate by effec- 
tively rejecting the /^-dependence of the real emission dijet contributions to the hadronic 
NLO W^y+jet cross sections. However, given that extra jet radiation becomes important in 
the tails of the transverse momentum distributions, Fig. 4, the veto introduces substantial 
uncertainties in this very phase space region. This can be inferred from Fig. 5, where we ex- 
emplarily examine the impact of the fixed-scale variation fx F = /j, R = 50 GeV, . . . , 200 GeV 
on the photon's transverse momentum for both inclusive and exclusive production. For 
completeness we note that dynamical scale choices, such as fiR = [ip — maxp^, result in 
quantitatively similar uncertainty bands. While the distribution's uncertainty band's rela- 
tive size is uniform over the entire range of the distribution for inclusive production, the jet 
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Table 1 : Inclusive next-to- leading order cross sections and total K factors for the processes pp — > 
e + ^ e 7j + X and pp — >■ e^u e ^fj + X at the LHC for identified renormalization and factorization 
scales, hr = hf = 100 GeV. The cuts are chosen as described in the text. 
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Figure 4: Transverse momentum of the W and maximum jet transverse momentum at leading order 
(dashed) and next-to- leading order (solid). The lower panels display the differential K factor, with 
the K factors of the total inclusive and exclusive production plotted as dashed horizontal lines. 
The W momentum is reconstructed from its decay products in the transverse plane. 
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Table 2: Exclusive next-todeading order cross sections and total K factors for the processes 
pp — > e + v e ^j + X and pp — > e~v e ^j + X at the LHC for identified renormalization and factorization 
scales, fiR = fiF = 100 GeV. The cuts are chosen as described in the text. 



veto stabilizes da^y/dp^ exclusively in threshold region, which dominates the integrated 
exclusive rate. In fact, the distributions for fi = 50 GeV and /i = 100 GeV, which are used 
to generate the uncertainty band in Fig. 5, intersect at ]P T ~ 100 GeV signaling an accidental 
cancellation of the renormalization and factorization scale dependence for the small total K 
factors of the exclusive sample. The large uncertainty in the exclusive p^, distribution's tail 
then translates into an only mild overall scale dependence 0^1° • Similar conclusions have 
been drawn for H/Z+jet production in Ref. [8]. 

This is yet another example of the well-known fact that total K factors and total scale 
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Figure 5: Transverse momentum of the photon for exclusive (left panel) and inclusive W~^f + jet 
production (right panel). The uncertainty bands refer to fixed scale variations around /zr = [if = 
100 GeV by a factor two in the NLO computation only. 



variations tend to be misleading when quantifying the impact of QCD quantum corrections 
to a given process. A better understanding of the QCD effects can be gained from differential 
K factors of (IR-safe) observables O, 

which encode the phase space dependence of the corrections, projected onto the respective 
observable. With unsuppressed extra jet radiation at the LHC's large center-of-mass energy^, 
the distributions' shapes are highly altered when including NLO inclusive corrections. While 
the shapes of purely electroweak distributions, i.e. distributions of observables that involve 
only photon, lepton and missing energy, survive to NLO QCD for the most part of the 
phase space, semi-hadronic observables get significantly modified with respect to their LO 
approximations due to additional jet radiation. Representatively, we show the azimuthal 
angle-pseudorapidity separation of lepton and photon and the minimal distance between 
jet and lepton in Fig. 6. We also plot the azimuthal angle between lepton and photon and 



^The proton is probed at small momentum fractions x ~ 0.02 at LO, which results in 
(T LO (VK ± 7j)/a LO (M^ ± 7jj) ~ 2 for the chosen selection criteria. Note that this is qualitatively different 
from the situation encountered in NLO 14^7 production. 
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Figure 6: Lepton-photon separation and minimum jet-photon separation at LO and NLO. The 
horizontal lines display the K factors of the total inclusive and exclusive production. 



W and photon in Fig. 8. The purely leptonic observables also receive sizable modifications 
at the edges of the phase space, which are determined by the chosen cuts. If, e.g., the 
Wj system recoils against additionally emitted partons at NLO, the W and the photon are 
forced to larger rapidity differences (Fig. 7), which communicates to the azimuthal angle- 
pseudorapidity separation of lepton and photon at large values. Obviously this effect cannot 
be buffered by the additional jet veto. It is important to note that the pseudorapidity 
difference of the W and the photon in Fig. 7 is experimentally not observable because the 
neutrino's longitudinal momentum explicitly enters the observable's definition. We show 
the distribution for comparison only. All other observables do depend only on the missing 
transverse momentum, which is experimentally reconstructed from the event's calorimeter 
entries [37]. 

The typical signatures of W^+jet production at the LHC are dominated by configurations 
close to the p^-thresholds; the entire event is central with lepton and photon preferably 
emitted at small angular distances in the transverse plane. The collinear photon-lepton 
singularity is cut away by the requirement of Eq. (15). The photon is typically emitted 
collinear to the W in the transverse plane. The jet recoils against the W^y pair, and is 
back-to-back to the W and the photon in the azimuthal angle distribution. For events with 
W and 7 back-to-back (e.g. \<pwy\ > 150 deg, where the inclusive NLO cross sections drops 



13 




0.5 




0.5 



2 

\yw - V-y\ 



Figure 7: Lepton-photon and VF-photon (pseudo) rapidity separation LO and NLO. The horizontal 
lines display the K factors of the total inclusive and exclusive production. The W four-momentum 
is defined from its decay products: pjf = p^+f ■ 



by about 80% to cr^ = 131 fb) the jets tend to balance each other at small rapidity gaps 
of order one. These gaps are due to the dominant qg- and ^-induced partonic subprocesses 
of the dijet contribution. Due to the relatively large available phase space for additional jet 
emission, the corrections are particularly large for this phase space region, K = 1.6. This 
is in accordance with the large differential K factor around cf>wy ~ ^ m Fig- 8. Note that 
this is also a part of the phase space which is more sensitive to the distinct photon isolation 
scales. In addition, with W and photon back-to-back, events with \4>wy\ > 150 deg can be 
understood as "genuine" pp — > W'j + X events, subject to anomalous couplings' studies". 

Anomalous couplings generically modify the p\ distribution at large values, Eq. (10). 
It is therefore worth commenting on the impact of the NLO corrections onto the region 
of phase space characterized by very large p^ already at this point. Considering energetic 
events in the tails of the px distributions (e.g. p\ > 1 TeV) the picture is quite different 
from the situation we have described above. Jet emission is logarithmically enhanced in the 

"These configurations give rise to large momentum transfers in the trilinear WWj coupling. The QCD 
corrections are therefore important to understand the deviations that result from anomalous couplings. We 
will discuss this in more detail in Sec. 3.4. 
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photon and lepton, and photon and W. The horizontal lines display the K factors of the total 
inclusive and exclusive production. 



dominant gluon- induced subprocesses qg — >■ W'yQ already at LO, which can easily be seen 
from the Altarelli-Parisi [38] approximation of collinear emission q — >■ QW as demonstrated 
in Ref. [11] 

2 7 2 

da ( qg WlQ ) = da ( qg ^ qi) \ i og 2 Pt_ ? (2Q) 

l&Ti sin W m^r 

for a diagonal CKM matrix. The preferred situation is therefore a collinear W-jet pair that 
recoils against the hard photon. This region of phase space receives sizable QCD corrections: 
The extra parton emission in these events has no preferred direction in the azimuthal angle 
and is kinematically unsuppressed. The pjf distribution receives sizable corrections for the 
same reason. The uncertainties in this region of phase space are dominated by the dijet 
contribution, which is only determined to LO approximation in our calculation. Hence, our 
NLO correction does not improve the cross sections' stability in this extreme region of phase 
space. However, as the anomalous couplings affect the p^, distribution at values much lower 
than A ~ TeV, see Eq. (10), the NLO corrections give rise to perturbatively predictive 
deviations from the SM, see Sec. 3.4. 

Vetoing the second resolved jet, Eq. (18), removes most of the characteristics of additional 
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Figure 9: Leading order and next-to- leading order max(p^) and p2, distributions at the Tevatron, 
including the respective differential K factors. The horizontal line represents the K factor of the 
total inclusive production. 



jet radiation. The total exclusive cross sections are tabulated in Tab. 2. Comparing the 
uncertainty band of the photon transverse momentum in Fig. 5 and the respective phase 
space dependence of the QCD corrections in Fig. 4, we conclude that perturbative stability 
against variations of factorization and renormalization scales of dijet-vetoed pp — > iwyj + X 
production shows up as a subtraction of a O{o i a 2 s ) leading-order contribution (the real 
emission dijet contribution) from a relatively stable inclusive NLO prediction. The vetoed 
contribution is kinematically well-accessible and unsuppressed by QCD dynamics. The larger 
scale dependence of the vetoed distributions' tails compared to inclusive production remains 
as an echo. At larger values of pj,, exclusive production does not yield a perturbatively 
reliable result, which is also indicated by negative weights. 

3.3 NLO QCD Wjj production in the SM at the Tevatron 

For Tevatron collisions we find a total cross section of 

a NLO = (14.86 ± 0.03) fb (# = 1.35). (21) 
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Tevatron, including the respective differential K factors. The horizontal line gives the K factor of 
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The proton and the antiproton are tested at x ~ 0.2 so that gQ-hiduced subprocesses domi- 
nate the total hadronic cross section. The considerably lower center-of-mass energy compared 
to the LHC in combination with the cuts on the jet transverse momentum effectively intro- 
duces a jet veto, so that the Tevatron shapes resemble the NLO exclusive LHC distributions. 
The shapes of the transverse momentum distributions at large pt are overestimated by the 
LO approximation, yielding differential K factors of order 0.5 at NLO for the distributions' 
tails. Since additional jet radiation is kinematically suppressed compared to the LHC, the 
semi-hadronic observables typically receive smaller relative corrections around the cx NLO /cr LO 
rescaled LO distributions. Yet, QCD-radiation effects are still sizable, and events tend to 
be re-distributed to smaller minimum separations of the hadronic jets with respect to the 
lepton and the photon, Fig. 10. 
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the LHC for anomalous input parameters 
|Ako|,|Ao| < 0.5 with dipole form factor 
n = 2 and the cutoff scale A = 1 TeV. 



3.4 NLO QCD W7+jet production with anomalous WW7 cou- 
plings at the LHC 

We now include anomalous couplings to the NLO Wj + jet cross section predictions. Most 
stringent bounds on anomalous WW'y couplings are currently given by the combined analysis 
of LEP data of Ref. [17], 

1 + Ak = 0.984l°;{$ , A = -0.016l°;°g , (22a) 

and recent fits at hadron colliders are from the Tevatron D0 experiment Ref. [18] 

1 + Akq = 1.0711$ A = -0.0l°;°H . (22b) 

Both bounds are at 68% confidence level, extracted from data assuming A = 2 TeV and 
dipole profiles. Note, that both experiments are consistent with the SM prediction Akq = 
Ao = 0. Generically, these bounds select a region in the parameter space where the QCD 
corrections are particularly important. This can be inferred from a scan over a wide (and 
experimentally ruled-out) range of anomalous parameters in Fig. 11. We choose cuts in 
resemblance to the selection criteria that are typically applied by the ATLAS collaboration 
to probe anomalous trilinear couplings, e.g., Ref. [10] 

pi > 100 GeV , p e T , f T > 25 GeV , R h > 1.0 , (23) 

where denotes the missing transverse momentum. In addition, we choose inclusive 
hadronic jet cuts 

p> T > 20 GeV , S = 0.4 , R jt > 0.2 . (24) 

These cuts yield a too low total rate at the Tevatron to be phenomenologically important. 
This can also be inferred from comparing the Pj, distributions at the LHC and the Tevatron 
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applying the cut of Eq. (3.15). We show the W _ 7+jet cross section for parameters Ko,Ao that 
are consistent with the LEP bounds. We choose dipole form factors n = 2 and a cutoff scale 
A = 2 TeV [10,18]. 



for large values, where the effects of anomalous couplings will be visible, Figs. 5 and 9. We 
therefore focus on anomalous couplings at the LHC. The qualitative reason why the QCD 
corrections turn out large for parameter choices in the vicinity of the SM is easily uncovered 
by examining the corrections' p\ dependence. From Fig. 4 we infer** that K{pl r ) > K in the 
threshold region and K{jP T ) < K in the tail of the distribution. Consequently, the region 
of phase space, where the anomalous couplings' impact is well-pronounced, i.e. p^ % A, 
provides a smaller fraction to the NLO cross section for inclusive cuts compared to the LO 
approximation rescaled by the total K factor of inclusive production. Additionally, at low 
transverse momenta, the distributions are dominated by SM physics due to small momentum 
transfers in Eq. (10), so that they are largely independent of Ak , Ao in this particular phase 
space region. In total, not only a large fraction of the cross section, but also a large share 
of its increase compared to LO, is insensitive to the underlying anomalous parameters for 
experimentally allowed values Ak , Ao- This is completely analogous to anomalous W Z+]et 
production [9]. The NLO inclusive cross section is therefore less sensitive to the anomalous 
couplings than the LO cross section, and K = a NLO /a LO is large in regions, where the 
distributions are dominated by their low px behavior: K peaks around the SM, An = A = 0, 
Fig. 11. From Fig. 5 and the related discussion, it is also apparent that, in addition to 
lost perturbative stability for large pj, the effects of anomalous couplings are suppressed in 
exclusive W^y+jet production. 

To increase the sensitivity to k and Ao in the experimentally allowed range, we addition- 
ally require the W and the photon to be back-to-back in the transverse plane, by imposing 
an azimuthal angle 

\(/>Wj\ > 150 deg. (25) 

**We slightly abuse the notation: K(p^) is the differential K factor in the sense of Eq. (19) and K without 
parentheses refers to the total K factor of exclusive production. 
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This cut effectively mimics "genuine" W"f events with additional hadronic activity. Given 
the hard pt requirements, this selection criterion can be replaced by a cut on \<f>i 7 \ or AR^ 
without qualitatively changing the phenomenology (see also Figs. 7 and 8). The resulting 
variation of the integrated + jet cross section for parameters (Ak , A ) in the range of 
Eq. (22) is of order 10%, Fig. 12. Comparing this variation to the uncertainty inherent to the 
SM expectation at the given order of perturbation theory, which, e.g., yields a ~ 60.6 fb for 
(J>r = A*f = 50 GeV, we see that the cross sections' increase due to the anomalous couplings 
is compatible with the SM NLO scale uncertainty, signaling a vanishing sensitivity of the 
total rate to A/t , X . 

This, however, does not hold for differential distributions at large momentum transfers, 
e.g. for the p T spectrum, which receives large anomalous couplings-induced modifications 
of the distribution's tail. The altered spectrum is well outside the SM-uncertainty band 
for larger values of (Ak ,X ), with a particular sensitivity to A . Remember that A dials 
the dimension six operator in Eq. (8), which is not present in the SM. The characteristic 
enhancement vanishes when the anomalous parameters approach their SM values, and the 
shape deviations become comparable to the distribution's uncertainty. The larger cross 
section at large p T compared to the SM translates into an increased cross section for the 
back-to-back configurations, which is also visible in the pseudorapidity differences at small 
separation, Fig. 14. The anomalous couplings' impact on this distribution is qualitatively 
different from the QCD corrections, which exhibit K{p7 r ) < K for large Pj,. Therefore, 
the NLO cross section at small rapidity differences is smaller than the NLO-normalized LO 
distributions suggests, Fig. 7. Yet, the NLO uncertainty from integrating over the small 
Pt configurations cover the anomalous couplings effect entirely, already by varying the scale 
within a small intervall, as indicated in Fig. 14. Given the residual anomalous couplings- 
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induced deviations of the p^, shape, a more inclusive measurement strategy that relies on fits 
to the inclusive ]P T spectrum or on multivariate analysis of the distributions can supplement 
traditional techniques and appears to be practicable. This is also motivated by the overall 
theoretical uncertainty of order 10% becoming comparable to the estimated experiments' 
systematics at the reported order of perturbation theory. 



4 Summary and Conclusions 

In this paper, we have calculated the differential NLO QCD corrections to W'j production in 
association with a hadronic jet, including full leptonic decays and all off-shell effects of the 
W. We have given details on the calculation's strategy and have also discussed the effects 
of anomalous WW'j couplings. The corrections are sizable and exhibit substantial phase 
space dependencies, and should be included in phenomenological analysis which employ 
these processes, either as signal or as background. In particular, we have found the exclusive 
production's perturbative stability to be accidental — additional jet vetoing does not amount 
to a experimental strategy which is under good theoretical control at the given order of 
perturbation theory. Qualitatively identical results have been shown to hold for WZ + jet 
production in Ref. [8], where the dominant QCD corrections are identical and only get tested 
at different scales. This strongly indicates that the additional jet veto does not provide a 
meaningful procedure for the entire class of massive diboson plus jet cross sections beyond 
theoretical contemplation. 

Even if theoretically less favored due to kinematical obstruction, inclusive 1^7+jet pro- 
duction (and hence inclusive production) exhibits potential sensitivity to anomalous 
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couplings via shape deviations of the p\ distribution. Integrated cross sections for our inclu- 
sively chosen leptonic cuts, however, entirely loose their sensitivity to modifications of the 
electroweak sector due to the low-pr QCD uncertainties. Providing the NLO corrections to 
pp — > £ ± + 7 + f + jet + X, we realistically asses the impact of anomalous couplings on the 
characteristic p\ distribution in Fig. 13. For larger anomalous couplings that are still com- 
patible with the combined LEP measurements, the distributions significantly deviate from 
their SM expectation and fall well outside the SM distribution's uncertainty. Comparing to 
anomalous WZ+jet production, we find more sizable deviations in the distributions' shapes 
in the allowed parameter range. Whether the observed sensitivity can be carried over to the 
experiment represents a challenging question, which is beyond the scope of this work. We 
leave a more thorough investigation of this direction to future work. 
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